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PREFACE 

If experimental work is to be accompli shed involving 
the use of a nuclear reactor, accurate knowledge concerning 
the density of neutrons at various locations within the re- 
actor is a necessity, whether the reactor be critical or 
subcr i t ical . 

The neutron density distribution in a given reactor 
can be estimated mathematically by the use of several tech- 
niques including transport theory, multigroup diffusion the- 
ory, and single group diffusion theory. Of these various 
mathematical theories, the single group diffusion is the 
only one amenable to hand calculations. Unfortunately, how- 
ever, diffusion theory is also the least accurate of the 
three named, although it does provide a good first approxi- 
mation to the actual neutron density distribution within a 
reactor . 

In order to accurately determine the neutron den- 
sity distribution within a reactor, experimental measure- 
ments must be conducted using one of the various methods 
available for such measurements. One of the most convenient 



methods is that of foil activation measurements 



Control rods made of materials which can absorb 



great numbers of neutrons are used to control the neutron 
distribution within reactors. In critical systems,, such 
control rods are, of course, necessary for safe and proper 
operation of the reactor. 

The University of Texas maintains a subcritical 
reactor which is used for teaching purposes and for research. 
The purpose of this thesis is to accurately determine the 
effect of a pure cadmium control rod on the neutron density 
distribution in this reactor. 

The neutron density distribution is determined by 
foil activation measurements both with and without the cad- 
mium control rod. Single group diffusion theory is used to 
predict the distribution mathematically, and the theoretical 
distributions are compared with the experimentally determined 
ones. 

The experimentally determined neutron density dis- 
tribution without the control rod is compared with the dis- 
tribution with the cadmium control rod in place in order to 
clearly show the effect of the rod on the neutron density. 

The author wishes to express his appreciation to 
Dr. R. N. Little for his advice and guidance during his 
supervision of this research and to Dr. J. A. Scanlan for 



his many helpful suggestions concerning 
thesis. Sincere gratitude is also due 
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summing assistance with the equipment u 
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CHAPTER 



I 



THEORY 



A . One Velocity Group Diffusion Theory Approximation 
without Control Rods"^ 



In the single velocity group diffusion theory it is 



as sumed 


that all 


neutrons 


in the sy 


s tern b 


eing studied have 


the s 


ame 


velocity, whi 


ch 


in the pre 


sent c 


as e 


is 2200 meters 


per s 


e co 


nd, the 


veloci 


ty 


of thermal 


neu tr 


ons 


If the 


neu tr 


on 


populati 


on is 


de s 


cribed in 


terms 


of 


the neutron 


dens i 


ty, 


n(r,t). 


where 














n 


(r,t) = 


number 


of 


neutrons 


of spe 


ed 


v per unit 








volume 


around space 


point 


r 


at time t. 


then 


the 


neutron 


flux, 


0( 


r , t ) , is e 


xpr e s s 


ed 


in terms of 


n(r , t 


) by the usual de 


fin 


i tion 














0(r, t) 




vn(r,t) 









where, in this case, all values are for thermal neutrons 
and v is constant and uniform throughout the reactor. 

A verbal statement of the neutron balance equation 
for a differential volume element, dV, includes a term to 

^R. V. Meghreblian and D. K. Holmes, Reactor Analysis 
(New York: McGraw-Hill Book Co 0 , I960), Chap. 5. 
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account for the net gain or loss through the boundaries of 
the volume element; and is 



-Rate of increase 
of the neutron 
density in dV 
L about r 



-number o 


f 




r number of — 




r ne t 


number of 


neutrons 


pro- 




neutrons ab- 




neu 


trons 


es cap- 


duced in 


dV 




sorbed in dV 




ing 


from 


dV 


L per unit 


t ime J 




L per unit time J 




L per 


unit 


time J 


tatement ; 


in mathematical form ; becomes 


the 


equation 



|^dV = S(r,t)dV - Z a 0(r,t)dV - £(r,t)dV 



where S(r,t) is called the source 
number of neutrons produced per u 
term &(r,t) is called the leakage 
let loss of neutrons from dV due 
boundari es . 

In diffusion theory the le 



function and denotes the 
nit time in dV, and the 
function and denotes the 
to diffusion through the 

akage term becomes 



£(r, t) = - DV 2 0 (r, t) 



and the neutron balance equation becomes; under steady state 
conditions for which n(r,t) is constant; 
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- DV 2 0 ( r ) + L 0 ( r ) = S(r). 
8 . 



In this equation D is called the diffusion coefficient and 
is defined to he 



D = 



'tr 



3 [ Z n - L (cos 6) ] 

total s N ' avg J 



where \^ r is the transport mean free path for the medium 
containing the volume element dV for which the equation was 
derived. Z a is the macroscopic adsorption cross section for 
this medium, and S(r) is the source term. As indicated above 



^total Z s (cos 0) aV g] 

where Z_ is the macroscopic scattering cross section for the 
s 

medium, (cos $) aV g = 2/3A and is the average scattering angle 

for neutrons involved in scattering collisions with nuclei of 

the medium. , _ is the sum of Z and Z • 

total s a 

The net neutron current through a given area lo- 
cated at r is given by the vector equation 



J ( r ) = -DV0 ( r ) . 
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The partial neutron currents through an area are defined by 
the equations 



«j (r) = i 0(r) - R j V0 ( r ) | co s 9 

+ 42 n 

and 

j_(r) = 0 ( r ) + | V0 ( r ) | co s 0 n 

where 0 n is the angle between the unit vector normal to the 
area and the unit vector parallel to the vector V0(r). The 
difference of the partial currents is equal to J(r). In 
other words ; 



J(r ) = j(r) - j (r) . 



The boundary conditions which must be used in obtaining a 
solution for the above diffusion theory neutron balance equa- 
tion are as follows: 

(1) 0(r) must be finite and nonnegative in all regions 
of the medium in which the diffusion equation ap- 
plies. 

(2) At any interface between two media with different 
diffusion coefficients, the partial currents must 



be continuous through the interface. 



If there is 
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no external source at an interface, then the neutron 
flux, 0(r), and the net current, J(r), must he con- 
tinuous at the interface. 

(3) The flux 0(r) must vanish at some finite distance 
outside the boundary or else the partial current 
from the area outside the boundary must be zero. 

In diffusion theory, the distance beyond the re- 
actor boundary at which 0(r) = 0 is called the ex- 
trapolation distance, d (= 0.7104 ) . 

Two assumptions were made in this treatment con- 
cerning the reactor assembly without the control rod. It 
was assumed that the air gap along the central axis of the 
core did not exist. This air gap actually represented only 
approximately 1.20 per cent of the volume of the core; so 
this assumption was not considered unreasonable. It was 
also assumed that a point source of neutrons was located at 
the center of the reactor. This point source is treated as 
a boundary condition in the following mathematical develop- 
ment. 

In the case of finite cylindrical geometry the 
diffusion equation can be written 

- DV 2 0(p,z) + Z a 0(p,z) = S(p,z) 
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where S(p,z) = vE f /( p,z) for a multiplying medium and is 
zero for a nonmultiplying medium. is the macroscopic 

fission cross section for the medium and v is the average 
number of neutrons resulting from a single fission event. 

If both a core and a reflector are involved, a 
diffusion equation must be written for both regions as follows 

D c v2 ^ C (c , ^ z ) + v Zf c /> c (p> z ) - Z a J c (p,z) = 0 

Dp>V 2 j^ R (p,z) - Z a *S R (p,z) = 0 

R 

The solutions to these simultaneous equations must satisfy 
all the boundary conditions previously stated. 

The above equations can also be written 

V 2 jZ$ c ( p , z ) + B 2 ^ p, z ) = 0 
and V 2 jz$ R (p,z) - K^ R (p,z) = 0 

where B 2 = [vZ. - £ ]/ D and k„ e Z /D-q. In neutron dif- 

C X 0 ac o Jt\ aj^ 

fusion theory a term called diffusion length, L, is defined 



L 2 




such that 
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Using this definition in the terms above gives 



.2 1 



and 



B 



2 _ 1 



R 

rvz. 



L 2 



R 



- 2-1 



c L Z a c 



k co - 1 



where k^, = vZ.p / Z a is the infinite multiplication constant 
c c 

for a medium composed of the same material as the reactor 



core . 



In cylindrical coordinate the V 2 operator is 



V 2 = i. 
P 



dp 



dz 



2 > 



so the above diffusion equations become 



and 



I 

p sp 

i 

P dp 



d »! c (o,b 



P f^ R (P;> z ) 



dp 



+ |^2 J^ c (p,z) + B 2 ^ c (p,z) = 0 



+ fp K {p ’ z) - K l*V p ' z) = 0 



both of which can be solved by the method of separation of 
variables . 

For the core it is assumed that 



jMp> z ) = r c(p)M z ) 
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and the resultant general solutions to the separated dif- 
fusion equation for the core are. if B 2 = k 2 + k 2 , 

^ 7 c p z 7 

Z (z) = A sin k z + C cos k z 
c ' z z z z 

R(p)=AJ(kp)+CYfkp) 
c VH/ p o v p M/ p O v P M/ 

where J o (k^p) an( ^- Y Q (kpp) are zero order Bessel functions 
of the first and second kinds, respectively. The function 
Y 0 (k p p), which approaches negative infinity as p approaches 
the value of zero, becomes more negative for small values 
of (kpp) than the function J Q (kpp) becomes positive. Thus 
since the boundary conditions stipulate that j^(p,z) can 
never be negative, C p must be taken to be zero. Since 
sin k z z becomes negative when z becomes negative, this term 
would violate the physical requirement that ^(p,z) be sym- 
metric about the z = 0 plane. Thus, A z must also be taken 
to be zero. 

Hence, the solution for the core region is 

J^ c (p,z) = A c J 0 (k p p) cos k z z . 



For the reflector it is also assumed that 
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V p ' z) = r r (p)z r ( z) 

and the resultant general solutions to the separated dif- 
fusion equation for the reflector are. if k 2 = k 2 - k 2 , 

R P z 

Z (z) = A„ sin k z + C cos k z 
R Z z z z 

Rr(p) = A p I 0 U p p) + C p K 0 ( Kp p) 

where I 0 ( K p p ) and K G (K p p) are zero order modified Bessel 
functions of the first and second kinds, respectively. Here 
again, the condition of symmetry of the neutron flux about 
z = 0 requires that A z be taken to be zero. 

Hence, the solution for the reflector region is 

j^ R (p,z) = f ArIq ( K pP ) + c rK 0 (k p p)] cos k z z . 

If p Q and z Q are defined to he the extrapolated 
radius and height, then the flux must he zero at the points 
where p = p Q or where z = z QO This boundary condition re- 
quires that ^(p,z) = 0 when p = p Q or when z = z Q . So 
K z z o = n / ^ since cos n/ 2 = 0. Thus k z = Jt/2z o . 

This boundary condition also requires that 
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[A R I o ( Kp p) + C r K o ( k 0 D ^ = 0 Wh6n 0 = °o‘ 



Thus 



A R 



C K o ( K o° o) 
Iq( k d°o) 



The equations for the flux in the reactor thus become, in 
summary , 



/> c (p,z) = A c J 0 (k n p) cos k z z 



^R ( P > z ) - C 



R 



K o ( K 0 p o ) 

K o( k p p) " I 0 ( Kpp c ) I o( |C p p ) 



COS 



2z, 



where z Q and p Q are the extrapolated height and radius, re- 
spectively . 

At the interface between the core and the re- 
flector, the boundary conditions require that the flux and 
the net current be continuous across the interface. Since 
the current is a vector quantity, this condition requires 
that the components of the current vectors be continuous. 
Thus, at the interface denoted by height Zj_ and radius p, 

^ c ( p > z i ) = p ’ z i) 

J c (p> z i) = J R (P; z i) 
z z 



and 
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which leads to the equations 



a c j o(* d p) cos k z z i = C R 



K 0 (V) - Io(^pP) 



T ^ 1 VP' 



COS 



2 z. 



and 



Dc^eJo^np) sln k z z i 



2 z. 



d r c r 



K (k p„) 

K o (k pP ) - I o (^ Po ) I o( K pp) 



s in 



2 z. 



Dividing the second of these equations by the first gives 
the equation 



k z tan k z z^ 



D 



R 



2z o D c 



tan 



2 z. 



) 



which can be solved by numerical methods to determine k z 
after the values of d r> D q , z ^ , and z Q are computed or meas- 
ured for the reactor being studied. 

Once the values for k z and k z are known, the cor- 
responding values of k^ and can be computed from the re- 

lations 



B 



2 

C 



kf + 



k R = 



K 



2 

Z 



and 
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since B 2 and k 2 depend on the physical constants of the re- 
actor and these constants are known. 

The constant A is related to the constant by 
C n 

the interface boundary condition that the flux must be con- 
tinuous at the interface denoted by radius p^ and the height 
z o Thus (p^^z) = ^ (p ,z) or, if the point where z = 0 is 
taken, 




which will give Cr as a constant multiplied by A Q . 

The constant A q can be evaluated by the boundary 
condition of continuity of net neutron current at the inter- 
face between the external neutron source and the core. If 
comparison of theory with experimentally determined flux 
distributions is desired, then the maximum value of the flux 

can be used for A . 

c 

The values of D c , Dr, B c , kp, k 2 , kr, Kp, k 2 , and 
Cr for the reactor used in this experiment are computed in 
Appendix A. 

If the neutron density distribution, n(p,z), 
rather than the flux distribution, ^(p,z), is desired, the 
the constants A c and Cr are simply divided by the thermal 
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neutron velocity (220^000 cm/sec) to obtain the equations 



n (p,z) = A'J (k p) cos k z 
c M ' ' c o v p^ ' z 



n 



R ( P .> z ) = C E 



( k _ p _ ) 

MV> - fjffl I o ( «»'>> 

o p o 



COS K Z Z 



where A' c = A c /v th and C£ = C R /v th . 



B . One Velocity Group Diffusion Theory 
Approximation with Control Rods^ 

Since the diffusion theory differential equations 
are the same for the core and reflector regions regardless 
of the presence of a control rod, the general solutions for 
these equations remain the same as for the case when no con- 
trol rod is present. The effect of the control rod is ac- 
counted for hy the "boundary conditions used to evaluate the 
constants. 

The general solutions for the diffusion equations 

are 



/> c (p,z) = [A c J 0 (k p p) + C c Y 0 (k p p) ]cos k z z 



^R ( P .> z ) = [ArIq ( K pP ) + C rK 0 ( K p p) ]CPS K Z Z 



2 

R. V. Meghrehlian and D. K. Holmes, Reactor Analy - 
sis (New York: McGraw-Hill Book Co., 1960), Chap. 11. 
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where B 2 = k 2 + k 2 and kS = k p - k 2 as before, and J~(k^p) 
cpz Rpz ; o x p M/ 

and Y Q (k p) are zero order Bessel functions of the first 
and second kinds, respectively, and I q (k p) and K Q ( k p ) 
are zero order modified Bessel functions of the first and 
second kinds, respectively. 

The boundary conditions which were applicable for 
the case with no control rod remain applicable, and for the 
case of a black control rod fully inserted along the axis 
of a cylindrical reactor, the additional boundary condition 
that f) Q = 0 at the point where p = b Q must be imposed. The 
radius b Q is defined by the relation 

b 0 = b - € 

where b is the radius of the control rod and € is the extra- 
polation distance into the rod. 

When a black control rod is present, the flux be- 
comes zero for p >0, and so the term Y Q (kpp) need not be 
eliminated from the equation for ^ c (p,z). The boundary 
condition that ^ c (b Q ,z) = 0 yields the relationship 



^■c^o ( ^ p^o ) ~ ^ c^o ( k p^o ) * 



This relation can be used to eliminate C c and yield the 



equation 



’c (p ' z) = A c 



J (k Id ) 

J (k p) - 7 ° 7 T ~ rT Y ( k p) 
o v o' Y n ( k p b 0 ) ° 0 



cos k z z , 



where A c< ,k p , and k z are yet to be determined. 

The boundary condition that ^(p,z) = 0 at the ex- 
trapolated radius, p Q , and height, z Q , yields, as for the 
case of no control rod, the relationships 



and 




A R 



'R 



Kp ( K pPo ) 

I o^ D P 0 ) 



where k p can be determined from the definition 
when z Q is known. 

In summary, the equations for the flux can now be 

written 



fi c (o,z) = a 



o p o ' 



cos k z z 



(>v>(p,z) = C 



R 



r ( \ K o (K P p o ) t ( \ 

K o U P p) " TUTI 1(3 K P P 
o' p M o' 



COS K Z Z 



where k z , b Q , and k are known. 
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The boundary condition for continuity of flux and 
net neutron current are used to determine the remaining con- 
s tan ts . 

At the interface denoted by height and radius 
the relations 



j^ c (p.< z i) = )^R(p; z i) 

and J c (p,z ± ) = J R ( p, Zi ) 

z z 

must be true. These relations again lead to the equation 

coskz. cos k z . 

z 1 _ z i 

■^c^z s i n ^z z i ^R k z s i n K z z i 
which can be written 

k„ tan k„z^ = — ^ n tan -ff . - z^ 
z z 1 D c 2 z Q 2 z 0 1 

and can be solved by numerical methods in order to determine 

k . After k is determined, k can be found from the rela- 
z z 9 p 

tion = kp + which has previously been defined. 

The relationship between A c and Cr is determined 
by the boundary condition that c = j^r at the interface 



denoted "by radius and height z. 



the relationship 



1 7 

This condition yields 



C 



R = 



A, 



^o ^ i ) 





Y o ( ^pPi ) 



if the z = 0 plane is used for the evaluation. 

In a similar manner, the boundary condition that 

J c = at the interface denoted by radius pj_ and the 

P P 

z = 0 plane can be used to determine A c since this would giv 
two simultaneous equations containing A c and Cj^ 0 In order 
to correlate theory with experiment, A c can alternatively be 
taken to be the flux observed along the midplane at the dis~ 
tance p for which 



Jq ( kpP ) 




= 1 . 



The values of all constants except A c and Cp? are 
the same as when no control rod was present, The values of 
A c and for the reactor used in this experiment are com- 
puted in Appendix A. 
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If the neutron density distribution, n(p,z), is 
desired rather than the flux distribution, j^(p,z), the con- 
stants A c and C R are divided by the thermal neutron velocity 
(220,000 cm/sec) to obtain the equations 



n c (p,z) = 



A T 



J 0 ( k P p) 



J o( k p b o) Y (k 

^ p o ) ° 0 



cos k z z 



n R (p } z ) = 



K o( K pp) 



k o<vo> t , ; 

I u o ) : ° ( » p) 

O v p M o 1 



COS K z Z 



where A' c = A c /v th and C£ = C R /v th . 



CHAPTER 



I I 



NEUTRON FLUX MEASUREMENT BY THE ACTIVATION METHOD 
A . General 

The activation method provides one of the most con- 
venient and most accurate techniques for measuring neutron 
fluXo This method is based on the fact that many elements 
become radioactive when exposed to a neutron flux. For sev- 
eral elements, this activation occurs only for a narrow 
range of neutron energies, By measurement of the induced 
activity, the neutron flux responsible for the activation 
can be determined. 

To measure a flux, a thin foil of the most suitable 
material is first irradiated for a certain time interval at 
the desired location. The foil is then removed from the 
neutron flux and the induced activity is measured in an ap- 
propriate counting device. The relationship between the 
measured activity and the neutron flux must include correc- 
tion terms for such things as counter efficiency and depres- 
sion of neutron flux by the foil being used to perform the 
measurement . 
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In order to be suitable for this type of measure- 
ment, the foil material should have a large neutron absorp- 
tion cross section for the energy range to be studied. Also 
the decay scheme of the radioactive isotope produced by the 
neutron activation should be simple enough for easy measure- 
ment and interpretation, and the half life for decay of this 
product isotope should be such that irradiation and count- 
ing periods are not excessively long. 

For thermal neutron detection, the most frequently 

3 

used elements which have strong resonance absorptions and 
yield induced radioactivity of suitable half lives and radia 
tions are In, Rh , Ag. Au, I, and Dy. 

B . Use of Dysprosium Foils for Thermal 
Neutron Detection 

The dimensions of the subcritical reactor used in 
this project made necessary the use of a Van de Graaff ac- 
celerator to produce the external source of neutrons neces- 
sary for a subcritical reactor (see Chapter III). Due to 
the operational characteristics of the accelerator used, it 
was possible to maintain a constant source strength for only 

3 C. F. Bonilla, Nuclear Engineering (New York: 
McGraw-Hill Book Co., Inc., 1957), p. 129. 
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short periods of time, and, hence, the foil detector ma- 
terial used had to be such that maximum induced activity 
could be reached in the shortest period of time feasible. 

Of the foil materials (in, Au, and Dy ) available, dysprosium, 
which has a large cross section for thermal neutron absorp- 
tion as well as a low cross section for epithermal neutron 
absorption, best met this requirement. However, dysprosium 
does have one disadvantage in that approximately twenty 
hours must be allowed for adequate decay of the induced ac- 
tivity prior to reusing a particular foil for another meas- 
urement . 

Preliminary calculations showed that, of the sev- 
eral naturally occurring isotopes of dysprosium which could 
result in radioactive products when irradiated by neutrons, 
only the isotope Dy 164 need be considered. Contributions 
to induced activity from neutron absorption by isotopes 
Dy 156 and Dy 158 are negligible due to low percentages of 
abundance and small neutron absorption cross sections as 
well as the long half life of decay of the radioactive prod- 
ucts Dy 157 and Dy 159 . The contribution to measurable ac- 
tivity by the transformation of Dy 165 into the radioactive 
isotope Dy 166 was determined by calculation to be less than 
2 X 10 -14 per cent. These calculations were confirmed by 
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experimental determination of the half life of decay for 
the activity induced in natural dysprosium by neutron ac- 
tivation. The observed half life was exactly that of Dy 165 . 

Thus, only the activity induced through neutron 
absorption by the isotope Dy 164 need be considered in con- 
nection with neutron flux measurements. The reactions in- 
volved are: 



Dy 165 + 7 



Dy 164 + n - Dy 



16 5m 



Ho 165 +|3 +7 



Dy 164 + n -*• Dy- 100 -* He 103 + p + y 



165 _ Ha 165 



The decay scheme for Dy 165 is shown in Figure II-l. 

The general differential equation for formation 
of radioactive product nuclei by the neutron absorption re- 
action is, in word form, 



'Rate of accumula- 




* Rate 


of produc-’ 




'Rate of loss 


tion of radio- 


= 


t ion 


of radio- 


- 


of radio- 


-active nuclei 




- ac tiv 


e nuclei 




-active nuclei^ 



In the case of Dy 164 two differential equations are neces- 
sary since Dy 164 can be transformed into both Dy 165 and 
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Me v 



0. 985 



0.555 

0.515 



0 . 356 
0.275 



0.095 

0 



Figure II-l. 



Decay Scheme for Dy 165 (from 
Tornau* ) 



*Rtidiger Tornau, Zeitschrift fur Physik , 159 , p. 

101, 1960. 
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Dyi 65m by the neutron absorption reaction. These equations 



are : 




(1) 5ra 
v ’ dT 


- ' X 5m D 5m< T > 


d Dr 

(2 ) dT 


(T) 

- ^/Vth + °- 976 X 5m D 5,.< T > 
- X 5 D 5< T > - 


The notation 


used in the above equations is as follows: 


d 4 


= number of Dy 164 nuclei present (assumed to be 
cons tant ) 


d 5 (t) = 


= number of Dy 165 nuclei present at time T, 


D 5m( T ) = 


= number of Dy 1G5m nuclei present at time T. 


T 


= irradiation time 


^th 


= thermal neutron flux during irradiation 


= 


= thermal neutron absorption cross section for 
activation of Dy 164 to Dy 165 . 


^4-5m 

a ny 


= thermal neutron absorption cross section for 
activation of Dy 164 to Dy 165m . 


a 5 " 6 

u ny 


= thermal neutron absorption cross section for 
activation of Dy 165 to Dy 166 . 
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X. = decay constant for Dy 165 . 

o 

X = decay constant for Dy 165m . 

5m 



Preliminary calculations showed that the effect of 
the activation of Dy 165 to Dy 166 is completely negligible, 
so the solutions to these two simultaneous differential equa- 
tions are found, by standard methods, to be 



4-5m J) I r i 

D 5» (T) * 7 X 5m L 1 - 



and 



D 5 (T) = B 



1 - exp (~X.,_T) 



+ C 



exp (-X. 5m T) - exp (-X^T) 



where B = 



ct 4_5 D. + 0.97 6 u 4-5m D 

!U_i /th 



and 



C = 



0.976 o-n~ 5m D 4 
^5m ” 



^th • 



In determining these solutions the boundary condition that 
both D (T) and D 5m (T) were zero at time T = 0 was used. 

If the foil is removed from the neutron flux at 
time T, the number of Dy lS5 and Dy l65m nuclei present at any 
time t thereafter is determined by the solutions to the dif- 



ferential equations 
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and 




= 0.976 ^5 m D 5m (t) - A. 5 D 5 (t) . 



In this case the boundary conditions to be used are that 
®5m U) = D 5m (T) and D^(t) = D^(T) at t = 0 ; where t is time 
after removal from the neutron flux. 

The solutions to these simultaneous differential 
equations are 



product of the decay constant and the number of nuclei pres- 
ent at that time, the following equations give the activity 
as a function of the time after removal of the detector foil 
from the neutron flux: 




and 







+ D 5 (T) exp (- X. 5 t) 



Since the activity at any time t is simply the 



A 5m( t ) = x 5m D 5ra( t ) 



A 5 ( t ) = ^"5^5 ( t ) 
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Thus, if the activity is counted for a specific 
time interval, t^ to t^, the number of radioactive disinte- 
grations, N, occurring during this time will be 



N = I A 5m ( t ) dt + I A 5 (t) dt 



dt + / V>5 (t > dt 



where N = total number of disintegrations occurring between 
times t^ and tg . Evaluation of this integral yields 



N = 



0.976 X 



5m 



D ( T ) - -7— v u cz 

. 5 U 5 ‘ 5m 



D C „(T) 



exp (- X 5 ti) - exp (- X g t2) 



1.976 X, 5 - X. 5m 



“ ^5m^ 



D (T) 
5m v ' 



exp (-^ 5m t 1 ) - exp (- ^ 5m t 2 ) 



This equation for N appears to be extremely cum- 
bersome and unwieldy, but after the constants are evaluated 
and the time intervals to be used for irradiation and count- 
ing are decided upon, the above equation reduces to N = 
constant X m i> ^ w here m is the weight of the detector foil 



in grams . 
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Trial runs and calculations using foils of various 
weights indicated that an irradiation time of five minutes 
( T = 0 to T = 5 ) and a counting interval of 10 minutes 
(t = 15 to t = 25) would give very accurate results. Evalu- 

Ca 

ation of the various constants needed for the above equa- 
tions was carried out and the results are tabulated in 
Table II-l. Using these constants, the equation for N re- 
duces to 



N = 0.6402 mjZ$ th 



where m = foi 
neutrons per 
trons per squ 
thermal neutr 



1 mass in grams and = neutron flux in 

cm 2 -min. After converting the flux into neu- 
are centimeter per second, the equation for 
on flux can be written as follows: 



neutrons 



th ' 



cm 



secy 



N 

38 . 412 m (grams ) 



C. 



Determination of Thermal Neutron Flux 
from Measured Activity 



iThe measured activity must be multiplied or di- 
vided by several correction factors in order to account for 
losses occurring in the counting process, according to 
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TABLE II-l 

CONSTANTS TO BE USED IN EQUATION FOR NUMBER OF 
DISINTEGRATIONS DURING COUNTING INTERVAL 



Cross Sections* a^ -3 = 800 X 10 -24 cm 2 

nj 



a 4-5m = 2000 X 10~ 24 cm 2 
ny 



Decay Constants = — 8-883 — = 0.004978 min""''' 

13 139.2 min 



Product of Decay 
Constants and 
Times 

~T = 5 min” 

tj_ = 15 min 
_tg = 25 min_ 



Values of 
Exponentials 



^5m 



0.695 
1.25 min 



= 0.5544 min 



x 5 t = 


0.02489 


\ 5m T = 2.7720 


\ 5 tl = 


0.07467 


yi = 8 ■ 3160 


X 5 t 2 = 


0.12445 


^5m^2 = 13 • 8600 


-S T 


= 0.97541 


e 0.062536 


e^ st l 


= 0.92805 


e^Si^l = 0.000245 


e " x 5 t 2 


= 0.88298 


e “ X 5 m t 2 _ o.OOOOOl 


e ' X 5 t l 


- e"^ 2 


= 0.04507 



g-^nitl _ e - x 5m t 2 = 0 . 000244 



*Hughes, Magurno, and Brussel, Neutron Cross Sections , 
Supplement No. 1 to BNL 525 , U. S. Government Printing Office, Washing- 
ton 25, D. C., January 1960. 



